library(sf)
package ‘sf’ was built under R version 3.6.2Linking to GEOS 3.7.2, GDAL 2.4.2, PROJ 5.2.0
library(tidyverse)
library(ggplot2)
loc_authority_with_extra_density_2019 <- read_csv("data/loc_authority_with_extra_density_2019.csv")
urban_rural_2019_with_ex_den <- read_csv("data/urban_rural_2019_with_ex_den.csv")
la <- st_read("https://opendata.arcgis.com/datasets/2698b6825e88404ab6272a5ff7f59ee6_0.geojson")
loc_authority_with_extra_density_2019
loc_authority_with_extra_density_2019 %>%
filter(feature_code == "S12000044")
loc_authority_with_extra_density_2019 <-
loc_authority_with_extra_density_2019 %>%
mutate(feature_code = recode(feature_code, S12000049 = "S12000046", S12000050 = "S12000044"
)
)
urban_rural_2019_with_ex_den <-
urban_rural_2019_with_ex_den %>%
mutate(feature_code = recode(feature_code, S12000049 = "S12000046", S12000050 = "S12000044"
)
)
la %>%
filter(lad18nm == "North Lanarkshire") #S12000046
Simple feature collection with 1 feature and 10 fields
geometry type: MULTIPOLYGON
dimension: XY
bbox: xmin: -4.19449 ymin: 55.73476 xmax: -3.711127 ymax: 56.03136
CRS: 4326
objectid lad18cd lad18nm lad18nmw bng_e bng_n long lat st_areashape
1 365 S12000044 North Lanarkshire 277984 665608 -3.9514 55.8681 472257663
st_lengthshape geometry
1 160705.9 MULTIPOLYGON (((-4.020121 5...
la_loc_authority_ex_den <- inner_join(loc_authority_with_extra_density_2019, la, by = c("feature_code" = "lad18cd"))
la_loc_authority_ex_den
la_urban_rural_2019 <- left_join(urban_rural_2019_with_ex_den, la, by = c("feature_code" = "lad18cd"))
la_urban_rural_2019
ggplot() +
geom_sf(data = la_loc_authority_ex_den, mapping = aes(geometry = geometry, fill = a_5_minute_walk_or_less), colour = "grey")

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkoc2YpCmxpYnJhcnkodGlkeXZlcnNlKQpsaWJyYXJ5KGdncGxvdDIpCgpgYGAKCmBgYHtyfQpsb2NfYXV0aG9yaXR5X3dpdGhfZXh0cmFfZGVuc2l0eV8yMDE5IDwtIHJlYWRfY3N2KCJkYXRhL2xvY19hdXRob3JpdHlfd2l0aF9leHRyYV9kZW5zaXR5XzIwMTkuY3N2IikKCnVyYmFuX3J1cmFsXzIwMTlfd2l0aF9leF9kZW4gPC0gcmVhZF9jc3YoImRhdGEvdXJiYW5fcnVyYWxfMjAxOV93aXRoX2V4X2Rlbi5jc3YiKSAKYGBgCgpgYGB7cn0KbGEgPC0gc3RfcmVhZCgiaHR0cHM6Ly9vcGVuZGF0YS5hcmNnaXMuY29tL2RhdGFzZXRzLzI2OThiNjgyNWU4ODQwNGFiNjI3MmE1ZmY3ZjU5ZWU2XzAuZ2VvanNvbiIpCgpgYGAKCmBgYHtyfQpsb2NfYXV0aG9yaXR5X3dpdGhfZXh0cmFfZGVuc2l0eV8yMDE5ICAgICAjUzEyMDAwMDQ5ICAgICNTMTIwMDAwNTAKYGBgCgoKYGBge3J9CmxvY19hdXRob3JpdHlfd2l0aF9leHRyYV9kZW5zaXR5XzIwMTkgJT4lIAogIG11dGF0ZShyZWNvZGUoZmVhdHVyZV9jb2RlLCApKQpgYGAKCgpgYGB7cn0KbG9jX2F1dGhvcml0eV93aXRoX2V4dHJhX2RlbnNpdHlfMjAxOSA8LQogIGxvY19hdXRob3JpdHlfd2l0aF9leHRyYV9kZW5zaXR5XzIwMTkgJT4lIAogIG11dGF0ZShmZWF0dXJlX2NvZGUgPSByZWNvZGUoZmVhdHVyZV9jb2RlLCAgUzEyMDAwMDQ5ID0gIlMxMjAwMDA0NiIsIFMxMjAwMDA1MCA9ICJTMTIwMDAwNDQiIAogICAgICAgICApCiAgKQoKdXJiYW5fcnVyYWxfMjAxOV93aXRoX2V4X2RlbiA8LQogIHVyYmFuX3J1cmFsXzIwMTlfd2l0aF9leF9kZW4gJT4lIAogIG11dGF0ZShmZWF0dXJlX2NvZGUgPSByZWNvZGUoZmVhdHVyZV9jb2RlLCAgUzEyMDAwMDQ5ID0gIlMxMjAwMDA0NiIsIFMxMjAwMDA1MCA9ICJTMTIwMDAwNDQiIAogICAgICAgICApCiAgKQpgYGAKCgpgYGB7cn0KbGEgJT4lIAogIGZpbHRlcihsYWQxOG5tID09ICJOb3J0aCBMYW5hcmtzaGlyZSIpICAjUzEyMDAwMDQ2IFMxMjAwMDA0NApgYGAKCgpgYGB7cn0KbGFfbG9jX2F1dGhvcml0eV9leF9kZW4gPC0gaW5uZXJfam9pbihsb2NfYXV0aG9yaXR5X3dpdGhfZXh0cmFfZGVuc2l0eV8yMDE5LCBsYSwgYnkgPSBjKCJmZWF0dXJlX2NvZGUiID0gImxhZDE4Y2QiKSkKCmxhX2xvY19hdXRob3JpdHlfZXhfZGVuCmBgYAoKYGBge3J9CmxhX3VyYmFuX3J1cmFsXzIwMTkgPC0gbGVmdF9qb2luKHVyYmFuX3J1cmFsXzIwMTlfd2l0aF9leF9kZW4sIGxhLCBieSA9IGMoImZlYXR1cmVfY29kZSIgPSAibGFkMThjZCIpKQoKbGFfdXJiYW5fcnVyYWxfMjAxOQpgYGAKYGBge3J9CiMgTmVlZCB0byBhZGp1c3QgdGhlIHpvb20sIGFuZ2xlIGFuZCBjb2xvci4gCmdncGxvdCgpICsKICBnZW9tX3NmKGRhdGEgPSBsYV9sb2NfYXV0aG9yaXR5X2V4X2RlbiwgbWFwcGluZyA9IGFlcyhnZW9tZXRyeSA9IGdlb21ldHJ5LCBmaWxsID0gYV81X21pbnV0ZV93YWxrX29yX2xlc3MpLCBjb2xvdXIgPSAiZ3JleSIpCmBgYAoKCg==